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Abstract 

The experimental results on the falling down of scaled factorial moments 
in azimuthal variable (f> is studied in some detail. It is shown that this phe- 
nomenon may be referred to the influence of transverse momentum conser- 
vation. The existing experimental data from DELPHI, EMU08, NA22 and 
UA1 are successfully explained. Various methods are proposed to partly 
eliminate this influence and rule out the 'falling down' of factorial moments. 
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Since the last decade the anomalous scaling of scaled factorial moments (SFM's) or 



tence of dynamical fluctuation or multifractal structure of multiparticle spectrum in high en- 
ergy collision processes. The corresponding experiments have been performed in various kind 
of collisions — from e + -e _ , hadron-hadron, hadron-nucleus to nucleus- nucleus. Plenty of new 
and interesting experimental results have been put in front of us. Various efforts || |B| have 
been taken to understand the physics behind these experimental findings. 

In this respect, it merits attention that, in almost all the cases, the logarithm of SFM's 
versus that of the cell width 5 rises as the decreasing of S (or as the increasing of partition 
number M of phase space, where M = A/S, A is the total phase space region in considera- 
tion), when M is not very big. There is, however, a noticeable exception, i.e. the logarithm 
of SFM's with azimuthal angle <f> as variable fall down as the increasing of InM in the first 
few points || |^, || ^J, cf. Fig.l. This prevents us from properly understanding the physics 
behind the scaling behaviour of SFM's found in experiments. 

The aim of this paper is to study this phenomenon in some detail. We will show 
that the falling down of SFM's in azimuthal angle may be referred to the influence of 
momentum conservation (MmCn) constraint in the collision processes. Various methods will 
be implemented to partly eliminate this influence and explore the inherent scaling property 



As is well known [Q, intermittency refers to the power-law behaviour of the scaled fac- 
torial moments F q as function of a cell width S of the phase space, in which the hadron 
multiplicity is measured. F q is defined by 



where M is the number of cells in horizontal average, n m is the multiplicity in mth J-cell in 
a high energy collision event, and ( ) denotes vertically averaging over many events. 

However, it should be taken into account that the total momenta of collision system 
must keep invariant in the space-time evolution process of the high energy collisions. This 
momentum conservation constraint may create some correlation between the multiplicity 
n rn of distinct phase space cells. But what is the influence of this correlation due to MmCn 
constraint upon the values of F q for different partition number Ml In order to attack the 
solution of this question more directly, Let us at first take the 1-dimensional (ID) momentum 
phase space as an example. 

Let Ap denotes the total momentum district in consideration. The momentum pi of 
each particle is restricted by the conservation condition y\ —1 pj = 0, where n is the total 
number of particle inside the region Ap. 

Firstly, it is clear that MmCn has no influence on the value of F q for M = 1, because 
F q for M — 1 depends only on the total number n of particles in the whole region Ap no 
matter how these particles are distributed, cf. Eq.( Q). More briefly, F q is 'blind' to the 
different distributions of p, therefore it takes the same value regardless of whether there is 
or isn't MmCn. 




Q| has been studied extensivelyH with the aim of exploring the possible exis- 



of SFM's. 
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When M — 2 the influence of MmCn comes in force, diminishing the value of F q . This 
can be seen as the following. When M = 2 the total region Ap is divided into two bins 
with positive and negative p respectively. Let P+ and P_ be the probabilities for a particle 
to lie in these two bins respectively. If there were neither dynamical fluctuation DF (i.e. if 
the probability did not fluctuate in itself) nor momentum conservation constraint, these two 
possibilities would be the same: 

P + = P_ (in case of neither DF nor MmCn). 

In an individual event, the statistical fluctuation (SF) will make the multiplicity n + and n_ 
of particles with positive and negative p fluctuate, 

n+ 7^ n_ (randomly due to SF). 

Assuming SF to be of Bernouli or Poisson type, it can be eliminated by means of factorial 
moments method[Q, resulting in a flat SFM versus M, i.e. 

p(M—2) _ p(M-i) ^ o j jjgy-kgj. jyp nor M m Cn). 

If there is dynamical fluctuation, the P + and P_ themselves will be unequal in an individual 
event, resulting in much larger fluctuations in n + and n_, so that F q will rise from M = 1 
to M — 2. On the other hand, The momentum conservation constraint 5^*=iPj = tends 
to separate particles to the positive and negative p regions, which reduces the difference 
between n+ and rt_ caused by DF and/or SF, making F q smaller than what would be 
expected without the restriction of MmCn. 

Thus it is predicted that F q should not rise so steeply as required by DF or even may 
fall down from the first point (M = 1) to the second point (M = 2) in contradiction to the 
flat SFM required by pure SF, because of the correlation effect between multiplicity n + and 
n_ from MmCn. 

In Fig. 2a, b is shown the results of Monte Carlo simulation of a random cascading model 
with and without MmCn. In this model, the total region Ap is divided subsequently with 
the subdividing ratio A = 2. The probability w; for a subdivision is taken to be 

Wi = — x^-i-^ -, i = l,2, (2) 

where a is the parameter of fluctuation strength, < a < 1; 7"j is a random number in the 
interval [-1,1]. After v steps, the probability in fcth subdivided bin (k = %i . . . , i v }) is 

/IN (l) (2) («/) ,„s 

The n particles are then assigned to the subdivided bins according to the Bernoulli distri- 
bution 

n' 

P{ ni ,n 2 , ...,n K ) = — — -p(l)'>(2r • ■ -p{K) n *, (4) 

71i!7l2! • ' ' flK. 

where K = 2 y is the total number of bins. Taking randomly the momenta of the rik particles 
in the fcth bin (fc = 1,2, .. . , K), a Monte Carlo sample of particle momenta is obtained. 
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In the calculation we have taken the model parameter a = 0.5, the cascading has been 
performed up to v = 7th step. 

In this model the dynamical fluctuation is taken into account up to the vth step of 
random cascading, i.e. up to the bin width as small as bp = Ap/2". When v is large, bp is 
very small, the neglecting of further DF will cause trivial error in the final results. 

^From the construction of model the resulting SFM's have the property of anomalous 
scaling 

Fq(M) oc M v , (5) 

so that the log-log plot of hiF q versus InAf is a straight line pointing upward for the point 
of M = 2*f], cf. Fig.2a. 

In the above model there is, of course, not any restriction from MmCn in the resulting 
sample, i.e. 2lLi Pi 0- We use a simple method to put in the constraint of MmCn. Let 



p. 

n 
i=i 

and shift the momentum of every particle in the event by an amount p c 

Pi =Pi~ Pc- (7) 

Then the new momenta p[ obey the conservation law, 53i=iPi = 0- This transform does 
effect the distribution of p, cf. Fig. 3, but not very much, and at the same time the fractal 
structure existing in the original sample retains in the new one. Using the cumulative variable 
as usualj[l), the distribution becomes flat again and the SFM's can be calculated. 

The results are shown in Fig. 2b. It can be seen in the figure that the first point, 
corresponding to M = 1 has the same height as that of Fig. 2a while the second point, 
M = 2, is lower than that of Fig. 2a due to MmCn, resulting in a falling down of lnF 2 from 
the first to the second point, confirming the prediction given above. When M increases 
further, lni^ increases again due to dynamical fluctuation. 

Thus we see that the falling down of F q at the first few points is due to putting in 
force of the MmCn effect. If we take the positive and negative p particles as particles of two 
"events", calculating SFM's in these "events" separately and then take the vertical average 
over them, the falling down phenomenon will disappear. It is really the case in our model, 
cf. Fig. 2c. 

On the other hand, if we make a cut and consider only the particles with \p\ < p cut , 
then the remaining particles will no longer be restricted so strongly by MmCn and the falling 
down phenomenon may graduately disappear. In Fig. 2c? is shown the result of our model 
with p cut = 0.8, a straight line is recovered. 

Let us now turn to the more realistic case of 3-dimensional (3D) phase space. It is no- 
ticeable that the momentum conservation constraint of sample is very different in transverse 
and longitudinal phase space. Generally, not all final fragments of high energy collisions are 

lr The 'chain effect' may appear in the plot of ln_F 9 versus InM if the continuously diminishing scale is 
taken in this kind of random cascading model |10|. 
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considered in the intermittency experiments. In particular, the leading particles, which take 
away in the average 50% of longitudinal momenta, are neglected. Therefore, the restriction 
of longitudinal momentum conservation on the particles analysed is quite weak. In fact, the 
experimental data of SFM's in longitudinal phase space (e.g. rapidity y or pseudorapidity 
rj) do not fall down in the plot of h\F q ~ In M(y), cf. Ref.||]. 

On the other hand, the constraint effect of MmCn in transverse plan is a little different 
between the samples from the different experiments. It depends on the distinct experimental 
conditions, i.e. the analysed phase space district cover both the central and fragmental 
regions, or only the central region, or just part of the central region. 

There is usually no cut in transverse momentum in the fixed target experiments and 
the e + e _ collisions, and the transverse momenta of leading particles are almost zero. So if 
the produced particles in both central and fragmental region are considered (e.g. DELPHI 
collaboration^), the transverse momentum conservation must be satisfied very well. 

In some intermittency experiments, only final state hadrons of central region are anal- 
ysed. For example, the NA22 collaboration put a cut of rapidity \y\ < 2 in longitudinal 
direction|| although there is no cut in the transverse directions. So the transverse momenta 
of final hardons in fragmental region have been neglected. However, since many experimental 
finding had shown that the final particles in central and fragmental regions may come respec- 
tively from several relatively independent sources (or fire balls) |l2|, the transverse MmCn 
condition is also kept approximately for the particles of central region alone. 

In UA1 experimental data||, only a part of the particles in the central region are 
analysed. A pseudorapidity cut of \r]\ < 1.5 is taken in addition to the beam pipe cut 
of pt > 0.15GeV in transverse direction. So the influence of transverse MmCn in UA1 
experimental data is not so strong, which will manifest itself by the weaker falling degree of 
SFM's in the plot of lnF q ~ lnM(</>). 

In order to reveal the influence of transverse MmCn upon the anomalous scaling be- 
haviour of SFM's in 3D phase space, we generalize the above ID cascading model with 
MmCn into 3D phase space. Constructing a 3D random cascading model in the space of 
< y < 1, < pt < 1, < ({> < 2ir, the results of Monte Carlo simulation are shown in the 
first column of Fig. 5. The 3D (y,pt, 4>) curve is a straight line as required by the cascading 
model, and the 1-, 2-dimensional SFM's get saturated, which is referred to the well known 
projection effect Jl3[. 

Now let us impose the transverse MmCn condition on the random cascading model. 
After all the vector momentum pi of the produced particles are determined in a simulated 
event, we shift the vector transverse momentum p ti of every particle by an amount p tc — 



1 sr^n -* 




(8) 



Then the transverse MmCn condition 



it 




(9) 
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is satisfied. 

The distribution function of rapidity y is, of course, not influenced by this transform 
(not shown in figure). The effect of the translation upon the distribution of p t and 4> are 
shown in Fig. 4. It can be seen from the figure that the <f> distribution is almost uninfluenced 
by the transform due to the axial symmetry. So the cumulative variable has to be used only 
for p t . 

The results of calculation with MmCn put in in this way are shown in the second column 
of Fig. 5. It can be seen clearly that the falling down effect does take place for variable <j>. It 
is also interesting to notice that the minimum of SFM in the lni^ versus InM plot occurs 
at M(4>) — 4, contradicting the ID MC results with minimum location at M((j>) = 2. This 
is as expected, because M(cf>) = 4 means that the positive and negative transverse momenta 
in both directions (denoted by pt y and pt z respectively) are divided into different bins, while 
in the case of smaller M(<j>), M(4>) — 2 for example, the positive and negative transverse 
momentum in only one direction (say pt y ) are separated. Therefore, in the case of M(<f>) = 2, 
F q is 'blind' to the different distributions in p tz and the effect of reduction of DF by transverse 
MmCn is weaker correspondingly. This is why F q (M = 1) > F q (M = 2) > F q (M = 4). 

If we consider the four quadrants as four different "events" and then take the vertical 
average, the falling down phenomenon will disappear, as can be seen in the result of model 
calculation shown in the third column of Fig. 5. 

If we impose a cut on p t to weaken the influence of MmCn, and consider only particles 
with p t < Pf Ui , then the falling down effect gets weakened too and disappear at about 
p c t ut = 0.7, cf. Fig.6. 

Thus we have successfully explained the experimental finding on the falling down of the 
first few points in the lni^(i^) versus lnAf((/>) plot, cf. Fig.l. The transverse momentum 
conservation is satisfied in case that all the produced particles in both central and fragmental 
regions ^ are considered, or approximately satisfied if only those in central region are 
considered. According to the above argument the minimum of lnF q versus lnM(</>) should 
occur at about M((f>) = 4. It does occurs at about M(<j>) = 5 in the experimental data 
in Fig. la, b. On the contrary, in the pp collider experiments of UAl collaboration, only 
a part of the produced particles in central region are analysed, and there is a transverse 
\pt \ > 0.15GeV cut due to the existence of beam pipe. So the effect of transverse MmCn is 
much weaker, cf. Fig.lc and Fig. 6c. 

Although there is a difference between these two cases, i.e. in our model of Fig. 6c we 
have imposed a cut at the upper side of pt, while in UAl experiment (Fig.lc), the cut is 
imposed at the lower side, the physics is the same, both are the relaxation of the restriction 
of transverse MmCn due to the particles carrying only a part of transverse momenta. It is 
noticeable that, besides the weakening of the falling down of SFM, this effect also makes its 
minimum move to a smaller value of M, locating at M < 4. This is the case both for the 
Monte Carlo model and for the experimental data, cf. Fig.lc and Fig. 6c. 

We conclude that in order to explore the physics behind the experimental phenomena on 
the scaling behaviour of scaled factorial moments (intermittency phenomena), the influence 
of transverse momentum conservation must be taken into account. From the constraint of 
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transverse momentum conservation, the falling-down of SFM with azimuthal angle variable 
in DELPHI, EMU08, NA22 and UAl can be explained successfully. 

Using the method of "quadrant analysis", i.e. taking the particles in different quadrant 
as different "event", or imposing a pt cut, this influence can be partly eliminated, especially, 
the falling down of the first few points can be ruled out. So these methods are helpful to 
reveal the inherent scaling behaviours of a fractal system. 
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Figure Captions 



Fig.l Scaled factorial moments with azimuthal angle as variable, showing the 'falling down' 
phenomenon. Data taken from Ref's.|| [?], 

Fig. 2 The scaling behaviours of SFM's in one-dimensional random cascading model: 

(a) , original cascading model without MmCn constraint; 

(b) , the cascading model with MmCn considered; 

(c) , the result of 'quadrant analysis' to the cascading model with MmCn; 

(d) , the result of the cascading model with MmCn after a cut of \p\ < 0.8 is imposed. 

Fig. 3 Distribution of p in ID cascading model with and without MmCn. 

Fig. 4 The influence of shift of transverse momenta upon the distribution function of pt and 
4> in 3D (y,p t ,<p) phase space. The solid curve denote the distribution function after 
MmCn is inputted, and the dashed ones denote that without MmCn. 

Fig. 5 The scaling behaviours of SFM's and their lower dimensional projection in 3D cas- 
cading model, where the first row is results of 3D (y — p t — (f) phase space, the second 
and third rows are that of 2D (y — <fi) and ID (<fi) projections respectively 

(a) , the original cascading model without MmCn (the first column); 

(b) , the cascading model with MmCn (the second column); 

(c) , the result of 'quadrant analysis' for the cascading model with MmCn 
(the third column). 

Fig. 6 The results of the 3D cascading model with MmCn after a cut of \p\ < p1 nt is imposed, 
(a)-(d) correspond to different p^ ut values respectively. 
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